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Abstract 

This  paper  uses  a  large  elementary  reaction  mechanism  to  study  the  homogeneous  chemistry  of  methane  and  natural-gas  mixed  with  air 
and  steam.  Temperatures  and  residence  times  are  chosen  to  represent  SOFC  operating  conditions,  including  within  feed  lines  that  may  be  at 
elevated  temperature  but  without  any  electrochemical  or  catalytic  interactions.  Mole  fractions  of  six  major  species  (CH4,  02,  H20,  H2,  CO, 
and  C02)  are  presented  as  contour  maps  as  functions  of  temperature,  residence  time,  and  initial  fuel  mixture  compositions.  In  addition,  deposit 
propensity  is  predicted  by  the  sum  of  mole  fractions  of  all  species  containing  five  or  more  carbon  atoms,  designated  as  C5+.  Comparison  with 
chemical  equilibrium  predictions  shows  that  the  homogeneous  reactions  are  far  from  equilibrium.  These  results  indicate  that  the  composition 
of  the  fuel  mixture  entering  the  active  SOFC  region  might  be  significantly  different  from  that  originally  entering  the  fuel  cell. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  principle  solid-oxide  fuel  cells  (SOFCs)  can  operate  di¬ 
rectly  with  hydrocarbon  fuels,  although  in  practice  measures 
must  be  taken  to  inhibit  carbon-deposit  formation  [1-11].  It 
is  often  desirable  to  premix  steam  or  air  with  the  hydrocar¬ 
bon  to  counteract  deposit  formation  and  possibly  to  assist 
with  thermal  management  [11-13].  Depending  on  system 
design,  inlet  fuel  mixtures  can  be  exposed  to  high  temper¬ 
atures  within  flow-distribution  networks  before  entering  the 
active  regions  of  the  fuel  cell.  Although  gas-phase  chemistry 
continues  within  the  cell,  the  homogeneous  chemistry  in  that 
region  may  be  strongly  influenced  by  electrochemistry  and 
catalytic  heterogeneous  chemistry  (both  of  which  are  outside 
the  scope  of  this  paper). 

This  paper  uses  a  detailed  reaction  mechanism  to  model 
the  purely  homogeneous  pyrolysis  and  oxidation  chemistry 
associated  with  initial  mixtures  of  methane  or  natural-gas 
with  steam  and  air.  It  considers  temperatures  ranging  from 
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600  to  900  °C  and  residence  times  ranging  from  0. 1  to  1000  s. 
We  distinguish  between  methane  and  natural-gas  to  explic¬ 
itly  account  for  expected  increased  chemical  reactivity  of  the 
higher  hydrocarbons  commonly  found  in  natural-gas.  In  our 
calculations  we  assume  that  the  fuel  is  free  from  the  im¬ 
purities  (e.g.,  sulphur).  Homogeneous  chemistry  plays  two 
potentially  important  roles.  First  it  acts  to  convert  the  initial 
hydrocarbon  fuel  into  H2,  CO,  H2O  and  other  compounds. 
Thus  the  “fuel”  that  enters  the  active  regions  of  the  fuel  cell 
may  not  be  the  initial  fuel  (e.g.,  natural-gas)  supplied  to  the 
system.  Second,  molecular- weight-growth  via  homogeneous 
chemistry  can  ultimately  lead  to  polyaromatic  hydrocarbon 
deposits  [14], 

Results  are  presented  in  the  form  of  contour  maps  show¬ 
ing  species  mole  fractions  as  functions  of  temperature,  res¬ 
idence  time,  and  initial  fuel  mixture  composition.  Deposit 
propensity  is  predicted  by  the  sum  of  mole  fractions  of  all 
species  containing  five  or  more  carbon  atoms,  designated  as 
C5+.  Increased  temperature  causes  increased  reaction.  Gen¬ 
erally  speaking,  natural-gas  is  always  more  reactive  than 
methane  alone.  Fuel-air  mixtures  are  significantly  more  reac¬ 
tive  than  fuel-steam  mixtures.  By  comparing  with  chemical- 
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equilibrium  predictions,  it  is  evident  that  kinetic  process  con¬ 
trols  the  homogeneous  chemistry  under  conditions  relevant 
to  SOFC  operation. 

2.  Model  description 

The  model  is  based  on  the  constant-temperature-and- 
pressure  species  continuity  equation  [15] 

dFk  _  «kWk 
dr  p 

where  the  independent  variable  is  the  time  t  and  the  dependent 
variables  are  the  mass  fractions  Tk-  The  molecular  weights 
are  and  the  mass  density  p  is  computed  using  a  perfect-gas 
equation  of  state.  The  species  molar  production  rates  a>k  are 
determined  from  the  elementary  reaction  mechanism.  As  dis¬ 
cussed  below,  the  mechanism  involves  about  350  species  and 
3450  reactions. 1  All  reactions  are  reversible,  with  the  reverse 
rate  coefficients  computed  from  the  temperature-dependent 
equilibrium  constants. 

Mathematically  the  model  is  an  ordinary-differential- 
equation  initial-value  problem,  which  is  solved  computa¬ 
tionally  using  Dvode  [16].  The  Chemkin  software  is  used 
to  evaluate  chemical  production  rates  from  the  reaction 
mechanism  [17].  Although  a  constant- temperature  constraint 
is  imposed  here,  one  could  easily  make  other  choices.  For 
example,  one  could  impose  a  temperature-time  profile, 
such  as  might  be  derived  from  knowledge  of  feed-line 
temperatures.  One  could  also  solve  the  thermal-energy 
equation  to  determine  temperature  variations  associated  with 
the  reaction  kinetics.  By  simply  specifying  a  temperature 
the  results  can  be  generalized  into  sets  of  contour  maps 
without  specific  regard  to  a  particular  fuel-cell  system 
design. 

The  reaction  mechanism  used  here  is  an  extension  of  one 
used  earlier  to  describe  butane  pyrolysis  [14].  It  incorpo¬ 
rates  oxidation  and  pyrolysis  kinetics  for  hydrocarbons  up 
to  Ce-  Additionally  it  includes  the  molecular-weight-growth 
reactions  for  polyaromatic  hydrocarbon  (PAH)  formation  via 
the  hydrogen-abstraction-acetylene-addition  mechanism  of 
Frenklach  and  Warnatz  [18].  This  earlier  mechanism  is  re¬ 
ferred  to  as  the  “Sheng-Dean”  mechanism. 

The  mechanism  in  the  current  work  is  modified  from  the 
Sheng-Dean  mechanism  as  follows: 

1 .  The  rate  coefficients  for  hydrogen-abstraction  reactions 
from  cyclopentadiene  were  updated  using  recent  results 
for  the  standard  heat  of  formation  of  the  cyclopentadi- 
enyl  radical  [19].  The  new  value  is  261.2  kJ  mol-1,  some 
21  kJ  mol-1  higher  than  the  previous  value. 

2.  A  substantial  number  of  pressure-dependent  reactions, 
including  dissociation  of  stable  hydrocarbons,  beta- 


1  The  mechanism  and  the  associated  thermodynamic  database  are  avail- 
able  upon  request. 


scission  of  hydrocarbon  and  oxygenated  radicals,  and 
radical  recombination  reactions  were  updated  [20]. 

3.  The  methane-dissociation  reaction  rate  (CH4  =±  H  + 
CH3)  was  lowered  by  a  factor  of  4. 

4.  The  rate  for  hydrogen  abstraction  from  ethane  by 

methyl  +  CH3  ^  C2H5  +  CH4)  was  updated  to 

the  value  reported  in  the  GRI-3.0  mechanism  [21]. 


3.  Mechanism  validation 

There  is  limited  data  published  for  methane  oxidation  or 
pyrolysis  in  the  SOFC  operating  temperature  regime.  These 
data  include  methane  pyrolysis  near  1000  K  [22,23]  and  sim¬ 
ulated  natural-gas  oxidation  between  1100  and  1200  K  [24]. 
As  discussed  in  the  following  sections,  the  current  mech¬ 
anism  is  better  able  to  represent  the  data  than  the  earlier 
Sheng-Dean  mechanism. 

3.1.  Methane  pyrolysis 

Back  and  coworkers  [22,23]  measured  the  production  of 
C2H6,  C2H2,  H2,  C2H4,  C3H6,  and  C3H4  at  very  low  levels 
of  methane  conversion  at  1038  K  and  0.58  atm.  An  interest¬ 
ing  feature  of  the  pyrolysis  data  is  the  increase  in  rate  with 
time  (Fig.  1).  This  behavior  is  especially  evident  for  ethane, 
where  the  mole  fraction  is  essentially  constant  between  500 
and  1500  s  before  undergoing  a  significant  increase  beyond 
1500  s.  Fig.  1  shows  that  the  Sheng-Dean  mechanism  agrees 
well  with  the  ethane  production  at  shorter  times.  However  the 
current  mechanism  represents  the  ethane  data  much  better  at 
longer  times.  Fig.  1  shows  that  the  current  mechanism  does 
a  reasonable  job  of  representing  the  transient  behavior  of  all 
measured  species. 

3.2.  Methane  oxidation 

El  Bakali  et  al.  [24]  used  steady-state  jet-stirred  reactor 
experiments  to  study  the  reaction  of  methane/ethane/oxygen 
mixtures  at  various  temperatures  and  atmospheric  pressure. 
The  initial  composition  was  0.73%  CH4,  0.073%  C2H6, 
1.14%  O2  with  the  balance  nitrogen.  The  relatively  small 
ethane  addition  was  used  to  simulate  natural-gas.  Mole  frac¬ 
tions  of  CH4,  C2H6,  C2H4,  C2H2  and  O2  were  measured  at 
a  residence  time  of  0.14  s.  Only  the  lowest  temperature  data 
(1 100-1200  K)  is  used  here,  because  it  is  closest  to  the  SOFC 
operating  range. 

Fig.  2  illustrates  a  comparison  between  the  model  predic¬ 
tions  and  measured  mole  fractions.  As  with  methane  pyrol¬ 
ysis,  the  current  mechanism  represents  the  data  reasonably 
well.  However  the  model  does  predict  a  somewhat  stronger 
temperature  dependence  than  is  observed.  These  compar¬ 
isons  suggest  that  the  current  mechanism  should  provide  rea¬ 
sonable  predictions  of  gas-phase  kinetics  under  SOFC  oper¬ 
ating  conditions. 
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Fig.  1.  Comparison  of  the  model  predictions  for  mole  fractions  of  C2H6, 
C2H4,  C2H2,  C3H6,  and  H2  to  various  CH4  pyrolysis  experiments  at  1038  K 
and  0.58  atm.  [22,23].  Experimental  data  (a)  refers  to  the  tabulated  data 
reported  in  [22],  (b)  refers  to  data  extracted  from  Fig.  4  of  the  same  work. 


Fig.  2.  Comparison  of  the  model  predictions  for  mole  fractions  of  CH4,  O2, 
C2H6,  C2H4,  and  C2H2  to  the  CH4  oxidation  experiments  of  El  Bakali  et  al. 
[24]  done  at  1  atm. 


4.  Results  and  discussion 

The  modeling  here  is  motivated  by  the  fact  that  hydro¬ 
carbon  fuels  are  often  mixed  with  steam  or  air  prior  to  in¬ 
troduction  into  the  SOFC.  This  strategy  can  improve  system 
performance,  including  the  inhibition  of  deposit  formation. 
The  results  in  this  paper  are  based  on  modeling  the  tran¬ 
sient  homogeneous  kinetics  of  methane  and  natural-gas  re¬ 
acting  with  air  and  steam.  The  simulations  are  isothermal  and 
isobaric,  with  temperatures  ranging  from  600  to  900  °C  and 
pressure  fixed  at  1  atm.  Results  are  reported  for  residence 
times  ranging  from  0.1  to  1000  s.  For  the  fuel-steam  mix¬ 
tures  the  initial  fuel  mole  fractions  range  between  0  and  1. 
For  the  fuel-air  mixtures  the  initial  fuel  mole  fractions  range 


between  0.2  and  1.  An  initial  fuel  mole  fraction  of  0.2  cor¬ 
responds  approximately  to  an  equivalence  ratio  of  0  =  2.1. 
Mixtures  with  methane  mole  fractions  lower  than  0.2  would 
be  sufficiently  close  to  combustion  conditions  that  the  fuel 
would  simply  burn  and  thus  be  of  little  practical  interest  for 
SOFC  operation.  Even  though  there  is  little  practical  value 
to  SOFC  fuels  that  are  primarily  steam,  results  for  the  entire 
range  of  initial  fuel-steam  mole  fractions  are  presented. 

Although  methane  is  the  dominant  constituent  of  natural- 
gas,  it  is  well  known  that  the  non-methane  constituents  of 
natural-gas  can  cause  significant  differences  in  the  reaction 
chemistry.  For  purposes  of  this  paper,  we  assign  natural-gas, 
designated  here  as  CH^,  the  composition  given  in  Table  1, 
which  is  representative  of  pipeline  natural-gas. 
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Table  1 

Composition  of  natural-gas  (CHJ)  used  in  the  calculations 


Species 

Molar  (%) 

ch4 

94.5 

c2h6 

2.5 

C3H8 

0.3 

H-C4H10 

0.05 

/50-C4H10 

0.05 

n2 

1.7 

co2 

0.9 

The  results  are  presented  as  contour  maps  showing  mole 
fractions  of  the  six  major  species  as  a  function  of  residence 
time  and  the  initial  fuel  mixture  composition  at  a  fixed  tem¬ 
perature.  The  propensity  to  form  polyaromatic  deposits  is 
characterized  by  the  sum  of  mole  fractions  of  all  product 
species  that  contain  five  or  more  carbon  atoms,  designated 
here  as  C5+ . 

In  studying  butane  pyrolysis,  Sheng  and  Dean  found 
that  the  onset  of  deposit  formation  could  be  correlated  with 
predicted  levels  of  C5+  [14].  Many  of  these  high-molecular- 
weight  species  are  cyclic  unsaturated  compounds  that  are 
precursors  of  polynuclear  aromatic  hydrocarbons  (PAH). 
Moreover,  forming  the  first  aromatic  ring  is  generally  con¬ 
sidered  to  be  the  rate-limiting  step  in  the  formation  of  PAHs. 
Flowing  butane  in  quartz  tubes,  Sheng  and  Dean  reported 
deposit  formation  at  a  nominal  temperature  of  700  °C  and  res¬ 
idence  time  of  5  s.  Under  these  conditions,  they  predicted  that 
the  mole  fraction  of  C5+  was  approximately  10-3,  suggesting 
this  value  as  an  approximate  threshold  for  deposit  formation. 


4.1.  Fuel-air  mixtures 

At  600  °C  with  methane  as  the  fuel,  the  top  two  panels 
of  Fig.  3  (CH4  and  O2)  show  that  below  a  residence  time 
of  1  s  very  little  reaction  occurs,  except  at  low  mole  frac¬ 
tions  of  methane.  Beyond  1  s  the  mole  fractions  of  both  CH4 
and  O2  start  to  decrease.  By  10  s  nearly  all  the  oxygen  is  con¬ 
sumed  over  the  entire  range  of  mixtures,  producing  H2O,  H2, 
CO  and  small  amounts  of  CO2  (panels  3-6  of  Fig.  3).  Note 
that  once  oxygen  is  completely  reacted,  the  contours  become 
nearly  flat  (i.e.,  independent  of  time).  Fig.  3  shows  that  the 
H2O  and  CO2  mole  fractions  increase  monotonically  as  the 
initial  mole  fraction  of  CH4  decreases  and  the  mixture  be¬ 
comes  more  oxidizing.  However  H2  and  CO  mole  fractions 
peak  at  intermediate  levels  of  initial  CH4  mole  fractions.  At 
lower  initial  CH4  mole  fractions,  H2  and  CO  are  oxidized  to 
H2O  and  CO2.  As  the  initial  methane  mole  fraction  decreases 
toward  0.2  the  production  rate  of  all  product  species  is  de¬ 
creased,  likely  caused  by  a  lower  radical-production  rate  at 
the  lower  methane  concentrations. 

The  bottom  panel  in  Fig.  3  shows  predictions  for  C5+ 
levels.  With  such  low  levels  predicted  here,  there  is  little 
likelihood  of  deposit  formation  from  gas-phase  reactions  for 
any  fuel-air  mixtures.  However,  it  is  interesting  to  note  the 
increased  tendency  toward  molecular-weight-growth  at  in- 


Fig.  3.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  methane-air  mixtures  at  600  °C. 
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termediate  initial  methane  mole  fractions.  The  formation  re¬ 
actions  are  slower  at  low  air  levels  due  to  the  intrinsic  sta¬ 
bility  of  methane.  The  radical  pool  is  insufficient  to  drive 
the  molecular-weight-growth  kinetics.  At  high  air  levels 
the  molecular-weight-growth  reactions  are  inhibited  by  the 
lower  concentrations  of  methane. 

Fig.  4  summarizes  the  natural-gas-air  predictions  at 
600  °C.  Generally  speaking  the  kinetics  are  slightly  faster 
than  those  for  methane.  In  the  case  of  natural-gas,  oxy¬ 
gen  is  depleted  faster  at  low  mole  fractions  of  CT.  This  is 
due  to  the  fact  that  the  higher  hydrocarbons  present  in  the 
natural-gas  are  more  reactive  than  methane.  The  product  dis¬ 
tribution  is  similar  to  the  methane  case,  but  generally  re¬ 
actions  start  slightly  earlier.  The  most  significant  difference 
is  the  predicted  production  of  C5+ ,  with  substantially  more 
growth  predicted  for  the  natural-gas  case.  Such  a  result  is 
anticipated  because  the  natural  gas  contains  more-reactive 
hydrocarbons.  Moreover,  the  natural-gas  has  a  “head  start” 
toward  molecular-weight-growth  due  to  the  initially  present 
higher  hydrocarbons.  Note  however  that  the  overall  level  of 
molecular-weight-growth  is  still  predicted  to  be  very  low. 
Assuming  a  threshold  of  C5+  >  0.001  for  deposit  forma¬ 
tion,  the  likelihood  for  PAH  formation  at  600  °C  is  still  quite 
low. 

Figs.  5  and  6  show  predictions  at  700  °C.  For  methane-air 
mixtures  (Fig.  5),  there  is  relatively  little  effect  of  the  in¬ 
creased  temperature  at  high  initial  CH4  mole  fractions.  This 
may  be  seen  by  comparing  the  O2  contours  with  the  600  °C 
case.  However  at  lower  initial  methane  mole  fractions  re¬ 
actions  begin  earlier  than  at  600  °C,  which  is  evident  from 
comparison  of  the  major-product  profiles.  As  at  600  °C,  once 
oxygen  is  completely  reacted,  the  profiles  of  the  major  prod¬ 
ucts  become  nearly  constant  in  time.  The  major  differences 
at  700  °C  as  compared  to  600  °C  are  as  follows: 

1 .  H2  peaks  at  approximately  5%  for  both  temperatures,  but 
at  700  °C  the  peak  occurs  earlier  in  time  and  the  high-H2 
region  is  larger. 

2.  At  the  higher  temperature,  CO  continues  to  increase  as 
the  methane  mole  fraction  decreases. 

3.  The  mole  fraction  of  C5+  is  larger  by  almost  two  orders 
of  magnitude. 

Fig.  6  shows  predictions  for  the  natural-gas-air  case  at 
700  °C.  As  at  600  °C,  the  results  show  natural-gas  to  be  more 
reactive  than  methane.  The  relative  differences  at  700  °C  be¬ 
tween  the  two  fuels  are  similar  to  those  predicted  at  600  °C. 
Relative  to  methane,  natural-gas  leads  to  a  substantial  in¬ 
crease  in  the  predicted  C5+  levels. 

The  predictions  at  800  °C  are  shown  in  Figs.  7  and  8.  The 
methane-air  results  (Fig.  7)  indicate  that  oxygen  is  depleted 
even  earlier  at  higher  temperature,  generally  within  a  few 
seconds  residence  time.  Although  the  species  distributions 
are  similar  to  the  700  °C  cases,  the  time  dependence  is  now 
significantly  altered.  This  is  most  evident  for  H2,  where  the 
mole  fraction  sharply  increases  at  long  times  and  low  lev¬ 
els  of  air  mole  fraction,  long  after  oxygen  has  completely 
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Fig.  4.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  natural-gas  (CHp-air  mixtures  at  600  °C. 
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Fig.  5.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  methane-air  mixtures  at  700  °C. 


Fig.  6.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  natural-gas  (CHj)-air  mixtures  at  700  °C. 
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Fig.  7.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  methane-air  mixtures  at  800  °C. 


Fig.  8.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  natural-gas  (CHp-air  mixtures  at  800  °C. 
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reacted.  The  increase  in  Tb  is  accompanied  by  a  decrease 
in  CH4.  These  changes  result  from  the  fact  that  the  system 
begins  to  behave  autocatalytically,  similar  to  the  pyrolysis 
results  shown  in  Fig.  1 .  As  illustrated  in  Fig.  7,  this  increased 
H2  production  rate  begins  around  100  s,  which  is  somewhat 
earlier  than  observed  in  the  pyrolysis  experiments  (Fig.  1). 
Flowever,  one  expects  this  time  shift  since  the  calculations 
are  at  slightly  higher  temperatures  and  pressures. 

At  800  °C  continued  increases  in  C5+  levels  are  predicted. 
The  mole  fraction  is  now  higher  by  about  an  order  of  magni¬ 
tude  than  predicted  at  700  °C  at  any  particular  residence  time 
or  mixture  composition.  As  seen  at  lower  temperatures.  Fig.  8 
shows  that  the  natural-gas  is  more  reactive  than  methane. 
With  natural-gas,  products  appear  at  slightly  earlier  times 
and  the  deposit  likelihood  increases. 

The  predictions  at  900  °C  (Figs.  9  and  10)  indicate  that 
substantial  gas-phase  chemistry  occurs  long  after  oxygen  is 
depleted.  The  concentration  of  Fb  increases  dramatically  at 
long  residence  times  and  low  air  mole  fractions  as  methane 
continues  to  react.  The  mole  fraction  of  CO  is  generally  much 
higher  than  that  of  CO2 .  The  higher  temperature  significantly 
increases  the  tendency  to  form  deposits.  The  predictions  with 
natural-gas  as  the  fuel  are  quite  similar  to  those  for  methane. 
The  temperature  is  now  sufficiently  high  that  the  reactivity 
difference  between  methane  and  the  higher  hydrocarbons  is 
not  so  pronounced. 

4.2.  Fuel— steam  mixtures 

The  fuel-steam  mixtures  show  relatively  little  homoge¬ 
neous  reactivity,  even  at  the  long  residence  times.  Even  at 
800  °C  for  the  methane  case,  the  top  two  panels  of  Fig.  11 
(CH4  and  H2O)  indicate  almost  no  reaction.  The  natural- 
gas  mixtures  are  more  reactive.  As  shown  in  Fig.  12  some 
Fb  is  produced,  especially  at  low  initial  levels  of  H2O.  The 
natural-gas  mixtures  are  also  more  likely  to  produce  de¬ 
posits  at  long  residence  times  over  a  large  range  of  mixture 
compositions. 

Figs.  13  and  14  show  that  reactivity  at  900  °C  increases 
markedly,  especially  at  long  residence  times.  Hb  is  the  dom¬ 
inant  product,  and  steam  becomes  a  reactant.  Note  that  the 
predicted  hydrogen  levels  are  comparable  to  those  predicted 
for  the  air  case  at  long  times  and  low  air  levels.  At  900  °C 
even  the  methane-steam  system  can  experience  substantial 
gas-phase  molecular-weight-growth.  Nevertheless,  even  at 
900  °C,  there  is  virtually  no  reaction  at  residence  times  below 
10  s  for  the  methane  mixtures.  However,  natural-gas  mixtures 
are  much  more  reactive,  with  substantial  H2  production  pre¬ 
dicted  at  residence  times  below  1  s.  As  expected  for  natural- 
gas,  the  likelihood  of  deposit  formation  is  generally  much 
larger  than  for  pure  methane. 

At  900  °C,  long  times,  and  low  levels  of  H2O  mole 
fraction  there  is  an  increase  in  H2  accompanied  by  a  de¬ 
crease  in  CH4.  Similar  behavior  was  observed  in  the  methane 
and  natural-gas-air  mixtures.  These  changes  are  caused  by 
the  system  beginning  to  behave  autocatalytically,  similar 
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Fig.  9.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  methane-air  mixtures  at  900  °C. 
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Fig.  10.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  natural-gas  (CHp-air  mixtures  at  900  °C. 
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Fig.  1 1 .  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  methane-steam  mixtures  at  800  °C. 


to  the  pyrolysis  results  in  Fig.  1.  Under  these  conditions, 
pyrolysis  causes  increased  production  of  H2  and  higher 
hydrocarbons. 

To  explore  the  role  of  H2O,  the  predicted  profiles  of  CH4 
and  H2  for  80%  natural-gas  and  20%  FUO  mixture  are  com¬ 
pared  to  those  where  FUO  is  replaced  with  N2.  Fig.  15  shows 
that  the  CH4  conversion  and  H2  production  is  virtually  iden- 
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Fig.  12.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  natural-gas  (CHp-steam  mixtures  at  800  °C. 

tical  in  the  two  cases,  suggesting  that  hbO  has  little  effect  on 
the  pyrolysis  reactions  under  these  conditions. 

At  900  °C,  H2O  is  responsible  for  producing  small  levels 
of  CO  and  CO2  as  seen  in  Figs.  13  and  14.  The  relatively 
slow  endothermic  reactions  of  FbO  with  CH3  and  H  radicals 
provides  a  source  of  OH  via  the  following  reactions: 

(2) 


Fig.  13.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  methane-steam  mixtures  at  900  °C. 

H  +  H20  ^  H2  +  OH  (3) 

The  OH  can  in  turn  participate  in  the  following  chain  reac¬ 
tions  to  produce  CO  and  CO2: 

CH3  +  OH  ^  CH3OH  (4) 

CH3OH  +  H  ^  CH2OH  +  H2 


CH3  +  H20  ^  CH4  +  OH 


(5) 
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Fig.  14.  Predicted  mole  fractions  as  a  function  of  residence  time  and  initial 
extent  of  mixing  for  natural-gas  (CHJ)— steam  mixtures  at  900  °C. 


CH3OH  +  CH3  —  CH2OH  +  CH4  (6) 

CH2OH  ^  CH20  +  H  (7) 

CH20  +  H  ^  HCO  +  H2  (8) 

HCO  ^H  +  CO  (9) 

C0  +  0H^C02+H  (10) 


Fig.  15.  Comparison  of  CFL)  and  Hi  histories  for  an  80-20  mixture  of 
natural-gas  (CHJ)  and  HiO  and  an  80-20  mixture  of  natural-gas  and  N2 
at  900  °C. 

The  radical  concentrations,  especially  CH3,  are  lower  for 
the  steam  case  than  for  the  air  case.  Thus  these  reactions 
become  important  only  at  a  higher  temperatures,  despite  the 
fact  that  H20  concentration  is  higher. 

4.3.  Kinetics  versus  equilibrium 

In  the  fuel-air  cases,  many  of  the  contour  plots  show  prod¬ 
uct  mole  fractions  that  are  independent  of  time  once  oxy¬ 
gen  is  depleted.  It  is  interesting  to  know  if  this  steady-state 
behavior  is  the  result  of  chemical  equilibrium  or  not.  Thus 
direct  comparisons  are  made  between  the  kinetics  predic¬ 
tions  and  the  local  chemical  equilibrium.  The  constrained 
temperature-pressure  equilibrium  conditions  are  computed 
using  Chemkin  [17]. 

Consider  cases  for  natural-gas-air  and  natural-gas-H20 
mixtures  with  the  natural-gas  mole  fractions  equal  to  0.5. 
Figs.  16  and  17  compare  kinetic  predictions  at  residence 
times  of  1000  s  and  at  various  temperatures  to  the  equilib¬ 
rium  predictions.  It  is  clear  that  these  systems  are  far  from 
equilibrated. 

Fig.  16  shows  that  for  the  natural-gas-air  mixture  the 
largest  deviations  are  at  the  highest  temperature  where  the 
kinetics  are  fastest.  The  gas-phase  kinetics  leads  to  over¬ 
oxidation  of  the  fuel,  with  much  higher  mole  fractions  of 
C02  and  H20  than  expected  at  equilibrium.  The  eventual 
approach  toward  equilibrium  is  substantially  inhibited  since 
the  gas-phase  kinetics  of  reactions  of  either  H20  or  CO?  with 
H  atoms  (e.g.,  reaction  (3)  and  reaction  (10)  in  the  reverse 
direction)  with  the  hydrocarbon  species  are  very  slow,  pri¬ 
marily  due  to  the  fact  that  these  reactions  are  endothermic 
with  correspondingly  high  activation  energies.  From  Fig.  17 
it  can  be  seen  that  for  natural-gas-H20  mixture,  there  are 
differences  between  equilibrium  and  kinetic  predictions  over 
the  entire  temperature  range.  The  gas-phase  kinetics  are  very 
slow,  leading  to  much  lower  mole  fractions  of  CO  and  FI2 
than  expected  at  equilibrium. 

In  addition  to  major  species,  the  equilibrium  C5+  levels 
are  also  predicted.  Graphite  is  not  included  in  the  equilibrium 
calculations  since  it  is  unlikely  to  be  formed  homogeneously 
on  any  reasonable  timescale.  Fig.  18  is  a  map  that  shows 
the  regions  where  the  equilibrium  C5+  is  above  and  below 
a  mole  fraction  of  10-3.  These  results  suggest  significant 
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Fig.  16.  Comparison  of  kinetic  predictions  for  a  50-50  mix  of  natural-gas 
(CHJ)  and  air  at  1000  s  to  the  equilibrium  calculations  as  a  function  of 
temperature. 

composition  and  temperature  dependencies.  At  600  °C,  the 
air  mole  fraction  would  need  to  be  at  least  0.17  to  avoid 
deposits.  As  the  temperature  increases,  more  air  is  needed  to 
avoid  deposits.  By  900  °C  nearly  55%  air  is  needed  to  avoid 
deposits  at  equilibrium. 

It  is  interesting  to  compare  the  equilibrium  predictions 
with  the  kinetic  predictions  of  C5+.  At  600  °C  (Fig.  4)  pre¬ 
dicted  C5+  values  from  the  kinetics  are  significantly  lower 
than  10-3  for  all  conditions,  whereas  equilibrium  indicates 
that  if  CH4  mole  fraction  is  above  0.83  the  amount  of  C5+  is 
greater  than  10-3. 

At  700  °C,  equilibrium  predicts  that  C5+  exceeds  1 0  3 
at  initial  natural-gas  mole  fractions  above  about  0.74.  The 


Temperature(°C) 

Fig.  17.  Comparison  of  kinetic  predictions  for  a  50-50  mix  of  natural-gas 
(CHJ)  and  steam  at  1000  s  to  the  equilibrium  calculations  as  a  function  of 
temperature. 


Fig.  18.  Predictions  of  the  natural-gas  (CHJ)  mole  fraction  (in  natural-gas- 
air  mixtures)  at  which  the  equilibrium  C5+  exceeds  10-3  as  a  function  of 
temperature. 

kinetic  predictions  at  700  °C  also  indicate  that  at  very  long 
residence  times  the  C5+  exceeds  10-3  for  natural-gas  mole 
fractions  greater  than  approximately  0.7.  However,  from  the 
shape  of  the  contours  one  expects  that  at  residence  times 
greater  than  1000  s  the  predicted  C5+  levels  exceed  10-3  for 
natural-gas  mole  fractions  well  below  0.7.  The  reactions  that 
drive  molecular-weight-growth  are  substantially  faster  than 
those  that  will  ultimately  bring  the  system  back  to  equilib¬ 
rium. 

Kinetic  control  is  even  more  evident  at  800  °C.  Equi¬ 
librium  predictions  show  that  C5+  levels  exceed  10-3  for 
natural-gas  mole  fractions  greater  than  about  0.57.  However 
it  is  evident  from  Fig.  8  (bottom  panel)  that  significant  C5+ 
is  present  at  much  lower  natural-gas  levels.  For  example,  at 
residence  times  near  1000  s  the  kinetic  predictions  indicate 
natural-gas  mole  fractions  as  low  as  0.3  produce  C5+  levels 
above  10-3.  This  is  clear  evidence  of  an  equilibrium  over¬ 
shoot. 

At  900  °C  Fig.  18  shows  that  at  equilibrium  natural-gas 
mole  fractions  greater  than  about  0.46  are  needed  for  C5+  to 
exceed  10-3.  The  kinetics  results  are  much  different,  as  can 
be  seen  from  Fig.  10.  At  1000  s  the  C5+  level  exceeds  10-3 
for  natural-gas  mole  fractions  as  low  as  0.2. 

Fig.  19  show  the  regions  where  C5+  levels  exceed  10-3  as 
a  function  of  temperature  and  initial  natural  gas  in  a  natural- 
gas-steam  mixture.  The  kinetics  results  at  600  and  700  °C 
(not  presented  in  the  paper)  indicated  very  little  reactivity  in 
the  system,  thus  the  amount  of  C5+  is  always  well  below  10-3. 


Temperature  (°C) 

Fig.  19.  Predictions  of  the  natural-gas  (CHJ)  mole  fraction  (in  a  mixture 
of  natural-gas  and  steam)  at  which  the  equilibrium  C5+  exceeds  10-3  as  a 
function  of  temperature. 


446 


G.K.  Gupta  et  al.  /  Journal  of  Power  Sources  156  (2006)  434-447 


However  equilibrium  predictions  indicated  that  the  amount 
of  C5+  is  greater  than  10-3  for  mole  fractions  of  natural-gas 
greater  than  approximately  0.96  at  600  °C  and  0.9  at  700  °C. 

Fig.  19  shows  that  at  800  °C  the  equilibrium  C5+  exceeds 
10-3  only  for  natural-gas  mole  fraction  greater  than  approxi¬ 
mately  0.82.  However  the  kinetics  predictions  in  Fig.  12  show 
that  at  longer  residence  times  the  C5+  levels  exceed  10-3  for 
all  natural-gas  mole  fractions  greater  than  approximately  0. 1 . 

Equilibrium  predictions  at  900  °C  show  that  C5+  levels 
exceed  10-3  for  natural-gas  mole  fractions  greater  than  about 
0.73.  However  it  is  evident  from  kinetics  predictions  shown 
in  Fig.  14  that  C5+  exceeds  10-3  at  residence  times  around 
100  s  for  natural-gas  mole  fractions  below  around  0.1. 

Generally  speaking,  under  SOFC  conditions  the  homoge¬ 
neous  chemistry  is  not  in  equilibrium.  Even  at  times  as  long 
as  1000  s  and  temperatures  as  high  as  900  °C,  both  major 
product  species  and  C5+  levels  are  well  out  of  equilibrium. 
In  natural-gas-air  mixtures  the  departure  from  equilibrium 
results  from  the  dramatic  decrease  in  conversion  rates  once 
oxygen  is  depleted.  In  natural-gas-steam  mixtures  the  de¬ 
parture  results  from  the  very  slow  kinetics,  resulting  in  slow 
conversion  to  equilibrium  products.  The  molecular-weight- 
growth  kinetics  appear  to  be  faster  than  the  slow  homoge¬ 
neous  “gasification”  reactions  for  both  natural-gas-air  and 
steam  mixtures.  Thus  for  residence  times  of  interest  the  ki¬ 
netics  predicts  higher  levels  of  C5+  than  would  be  present  at 
equilibrium. 

5.  Summary  and  conclusions 

This  paper  reports  a  modeling  investigation  of  the  gas- 
phase  reactions  of  methane  and  natural-gas  mixed  with  either 
air  or  steam  at  temperatures  from  600  to  900  °C  and  at  a  wide 
range  of  initial  compositions  and  residence  times  relevant 
to  SOFC  operation.  The  results,  presented  as  contour  maps, 
provide  a  great  deal  of  information  summarizing  the  behavior 
of  homogeneous  chemistry  as  could  occur  in  non-catalytic 
feed  lines  even  before  fuels  enter  the  electrochemically  active 
regions  of  the  SOFC. 

Adding  even  small  amounts  of  air  to  methane  or  natural- 
gas  significantly  increases  conversion  to  H2,  CO,  fFO,  and 
C’Ot  .  Under  most  circumstances  oxygen  is  completely  con¬ 
sumed  within  10  s  at  600  °C  and  within  1  s  at  900  °C.  Thus 
the  fuel  composition  entering  the  active  SOFC  regions  can 
be  significantly  different  from  the  initial  room-temperature 
mixture  fed  to  the  system.  Such  conversion  can  significantly 
affect  the  cell’s  electrochemical  performance.  For  example, 
the  steam  that  is  formed  homogeneously  in  the  feed  lines  can 
participate  immediately  in  reforming  reactions  on  the  SOFC 
anode,  before  any  further  fFO  is  produced  by  charge-transfer 
chemistry. 

Steam  addition  to  the  fuel  mixture  causes  much  less  reac¬ 
tion  than  air  addition.  However  at  900  °C  and  residence  times 
above  approximately  100  s  there  is  significant  conversion  of 
the  fuel  to  H2  and  higher  hydrocarbons. 


In  addition  to  major  gas-phase  products,  the  model  also 
predicts  molecular-weight-growth  to  higher  hydrocarbons. 
The  sum  of  mole  fractions  of  all  species  containing  five 
or  more  carbon  atoms,  designated  C5+,  is  used  as  an  indi¬ 
cator  of  deposit  propensity  via  polyaromatic  hydrocarbons. 
The  methane-air  mixtures  generally  produce  much  higher 
C5+  levels  than  the  methane-steam  mixtures  at  lower  tem¬ 
peratures.  However  the  C5+  levels  in  natural-gas-air  mix¬ 
tures  is  comparable  to  that  in  natural-gas-steam  mixtures 
at  higher  temperatures.  At  900  °C,  C5+  mole  fractions  can 
exceed  10-3  at  residence  times  greater  that  a  few  sec¬ 
onds.  These  C5+  levels  raise  a  concern  about  deposit  for¬ 
mation.  Interestingly,  the  kinetic  model  predicts  C5+  levels 
well  above  equilibrium  levels  for  some  residence  times  of 
interest. 

As  should  be  anticipated,  higher  temperatures  lead  to 
faster  kinetics  and  natural-gas  reacts  more  rapidly  than 
methane.  Nevertheless,  the  qualitative  behaviors  are  similar. 
Some  of  the  largest  differences  between  methane  and  natural- 
gas  are  in  the  C5+  predictions,  with  natural-gas  being  more 
likely  to  form  PAH  deposits. 

There  are  regions  at  longer  residence  times  where  the  ki¬ 
netics  appears  to  be  in  steady  state  in  the  sense  that  species 
composition  changes  very  slowly.  However,  by  direct  com¬ 
parison  with  equilibrium  predictions,  the  kinetics  are  found 
to  be  far  from  equilibrated.  For  the  fuel-air  cases  the  system 
appears  to  be  kinetically  “over-oxidized”  in  the  sense  that 
products  like  H2O  and  CO2  are  greater  than  equilibrium  val¬ 
ues  while  H2  and  CO  are  lower.  In  the  fuel-steam  cases  where 
the  chemistry  is  very  slow,  the  product  species  are  predicted 
to  be  much  lower  the  equilibrium  values. 
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